function apso 
Lb=2.592e-3*[-7.8*5 -1.6*3 -3.8*6 -7.8*5 -1.6*3 -3.8*6 -7.8*5 -1.6*3 3.8*6 -7.8*5 -1.6*3 -3.8*6 -7.8*5 -1.6*3 -3.8*6 -7.8*5 -1.6*3 -3.8*6 -7.8*5 -1.6*3 -3.8*6 -7.8*5 -1.6*3 -3.8*6 -7.8*5 -1.6*3 -3.8*6 -7.8*5 -1.6*3 3.8*6 -7.8*5 -1.6*3 -3.8*6 -7.8*5 -1.6*3 -3.8*6 ];
Ub=2.592e-3*[4.8*5 7.2*3 3.2*6 4.8*5 7.2*3 3.2*6 4.8*5 7.2*3 3.2*6 4.8*5 7.2*3 3.2*6 4.8*5 7.2*3 3.2*6 4.8*5 7.2*3 3.2*6 4.8*5 7.2*3 3.2*6 4.8*5 7.2*3 3.2*6 4.8*5 7.2*3 3.2*6 4.8*5 7.2*3 3.2*6 4.8*5 7.2*3 3.2*6 4.8*5 7.2*3 3.2*6];
para=[1000 150 0.95]; 
[gbest,fmin]=pso_mincon(@cost,@constraint,Lb,Ub,para); 
 Bestsolution=gbest  
 fmin 